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An indirect adaptive control algorithm has been studied 
to control physical systems with parameter uncertainties. 
Although the particular algorithm investigated is applicable 
to a wide class of discrete time linear systems, parameter 
convergence and therefore global stability of the whole system 
is guaranteed only if the external excitation contains a 
sufficiently large number of frequency components. 

This research study has also investigated the possibility 
of stopping the identification procedure when the parameter 
error becomes sufficiently small, so the controller auto- 
matically turns itself from adaptive to time invariant, while 
still guaranteeing global stability of the closed-loop system. 
In this way the adaptive controller might be activated only 
when its gains do not provide satisfactory performances. 

Computer programs of the indirect adaptive control have 
been written for simulation purposes using both recursive 
least squares and projection algorithms. 
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I . 



INTRODUCTION 



One of the most challenging, interesting and active fields 
of Automatic Control is Adaptive Control. To implement high- 
performance control systems when the plant dynamics are poorly 
known or when large and unpredictable variations occur, the 
control engineers prefer to use an important class of control 
systems called Adaptive Control systems. Adaptive Control 
comes from a desire and need for improving performance of 
complex engineering systems with large uncertainties. It is 
especially important in systems with many unknown parameters 
that are changing with time. Also, it can be defined as a 
special type of nonlinear feedback control, as a nonlinear, 
nonautonomous dynamic system. 

An adaptive controller can change its behavior in response 
to changes in the dynamics of the plant and the disturbances. 

The term adaptive control has been used at least from the 
beginning of the 1950's. With recent advances in microproc- 
essor technology, it has become feasible to implement adaptive 
algorithms efficiently in real time at a reasonable cost. 

There are three schemes for parameter adaptive control : 
gain scheduling, model reference control and self-tuning 
regulators. The starting point is an ordinary feedback con- 
trol loop with a process and regulator with adjustable 
parameters. But the main problem is to find a convenient way 
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of changing the regulator parameters in response to changes 
in. process and disturbance dynamics. The following schemes 
differ only in the way the parameters of the regulator are 
adjusted . 

A. GAIN SCHEDULING 

Gain scheduling depends on finding auxiliary process 
variables correlated with the changes in plant dynamics. In 
this way it is possible to reduce the effects of parameter 
variations by changing the parameters of the regulator as 
functions of the auxiliary variables. 

The block diagram of this scheme is given in Figure 1.1. 




Figure 1.1. Block Diagram of Gain-Scheduling System 
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B. MODEL-REFERENCE ADAPTIVE SYSTEMS 



In this type of adaptive system, a reference model speci- 
fies the desired performance, and tells how the process output 
should respond to the command signal. A block diagram of 
model reference system is given in Figure 1.2. As seen, from 
this figure, the reference model is part of the control system. 




Figure 1.2. Block Diagram of the Model Reference 
Adaptive System 
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This adaptive system has two loops : the lower loop con- 

tains Regulator and Process. The upper loop adjusts the 
parameters of the regulator, in such a way that the error 
between the model output and the process output tends to zero. 
In this type of control system, the main problem is to deter- 
mine the adjustment mechanism so that a stable system results, 
which brings the tracking error to zero. 

Model reference adaptive systems can be further subdivided 
into two categories: direct and indirect. In indirect con- 

trol the plant parameters are estimated and the control param- 
eters are adjusted based on these estimates so that the overall 
plant transfer function matches that of the reference model. 

In direct control no effort is made to identify the plant 
parameters but the control parameters are directly adjusted to 
minimize the error between plant and model outputs. 

It turns out that the model reference approach is applica- 
ble only to plants with stable zeroes. In fact the only way 
the closed loop transfer function (Plant & Regulator) can 
match the one of the model in its poles and zeroes, is by re- 
moving the plant zeroes by cancellation with closed loop poles. 
This operation leads to the presence of uncontrollable or 
unobservable modes in the closed loop systems, which can be 
accepted only if they are stable (i.e., their effect decays 
to zero with time) , and this constitutes the major limitation 
for model reference adaptive systems. 
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C. SELF-TUNING REGULATORS AND POLE PLACEMENT 



A third approach to the adaptive problem is the self- 
tuning regulator. A block diagram of this type control system 
is given in Figure 1.3. 



Desired performance 




Figure 1.3. Block Diagram of a Self-Tuning Regulator 

The reference model of the previous approach is replaced 
by some more general desired performances; such as error 
minimization or desired closed loop poles. 

There are two loops in the system configuration. The 
lower loop contains the plant and linear feedback regulator. 
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The upper loop consists of a recursive parameter estimator 
and a design calculation adjusts the parameters of the 
regulator . 

Also it is possible to classify this adaptive control 
scheme as direct and indirect. This classification depends 
on the complexity of the design calculation block that is 
seen in Figure 1.4. 




Figure 1.4. General Block Diagram of an Adaptive 
Control System 
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The first method is to parameterize the system directly 
in terms of compensator parameters. 

This research project will concentrate on the last method, 
that is, pole placement and indirect adaptive control. The 
evaluation of the control law is indirectly determined on 
the basis of parameter estimates. This method is also called 
explicit, since the design is based on an explicit estimation 
of the process model.. The parameters of the controller are 
updated indirectly via estimation of the process parameters. 
Several algorithms to estimate the parameters of a linear 
model are available in the literature. In this thesis two 
well-known estimation algorithms will be investigated: Recur- 

sive Least Squares and Projection. They represent a tradeoff 
between complexity of computation and performances, in the 
sense that best performances (with relatively high complexity 
are obtained by using Recursive Least Squares. Also, Block- 
processing is investigated to determine the period of the 
adaptation of the control parameters. Finite time persistency 
of excitation is also studied and simulated for indirect 
adaptive control . 

The main reason of choice of indirect adaptive control 
for this research comes from the fact that it is applicable 
to nonminimum-phase systems, and there are no limitations on 
zeroes of the plant. Therefore, it is more general than model 
reference adaptive control. 

Because of the effects of the zeroes of sampled data 
systems on adaptive control theory, we start investigating 
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the behavior of zeroes of sampled data systems. This thesis 
is organized as follows: In Chapter II zeroes of sampled 

data systems are analyzed by stating results available in 
the literature, and simulation studies. The indirect adaptive 
control problem is presented in Chapter III and computer 
simulation studies are given in Chapter IV. 

The computer programs are presented in the Appendices B, 

C and D. Also the description of the programs is given in 
Appendix A. 
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II . ZEROES OF SAMPLED DATA SYSTEMS 

Many control strategies are based on the assumption that 
the zeroes of the plant are on a stable region of the com- 
plex plane so that they can be cancelled by precompensation 
or closed loop poles. One example is the model reference 
adaptive control mentioned in the Introduction. However it 
turns out that in sampled data systems with Zero Order Hold 
(the most popular mean of Digital to Analog conversion) some 
zeroes of the resulting discrete time plant are often in the 
unstable region. In this chapter we analyze the position of 
the zeroes of sampled data systems, in relation with the con- 
tinuous time plant dynamics and sampling frequency. The 
structure, the number of zeroes outside the unit disc and the 
low frequency characteristic of the pulse transfer function 
are examined from the transfer function for an important set 
of continuous time processes. 

Finally, it is investigated that zeroes of the sampled 
data systems are sensitive to high frequency poles present 
in continuous time transfer function. 

A. TIME DOMAIN ANALYSIS OF THE ZEROES OF SAMPLED DATA 

TRANSFER FUNCTIONS 

Poles and zeroes are important parameters of linear time- 
invariant systems. The zeroes describe the way the internal 
variables are coupled to the inputs and the outputs. As 
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described in [Ref. 6], the unstable zeroes limit the perfor- 
mance of control systems, since many design techniques are 
based on the cancellation of the process zeroes. This can be 
done provided the system's zeroes are stable. 

The transformation of poles in the continuous time domain 
to the corresponding discrete time is given as : 

P T 

P.. - e i (2.1) 

l 

where T is the sampling period. This transformation maps 
the left-half part of the s-plane onto the unit disc, so that 
stability is preserved. But there is no simple transformation 
for. zeroes from continuous to discrete time domain. The type 
of hold circuit affects the position of the zeroes. Most 
digital control systems use a zero-order hold, and for this 
type of hold circuit, the effects are considered. 

In the following discussion, the main results are limit 
theorems, which give the zero locations for small and large 
sampling periods. 

If the continuous time transfer function G(S) is rational 
it follows from Equation (2.2) that H(z) is also a rational 
function 

H ( z ) = (1-z -1 ) J ( (e P i T )/(z-e P i T ) )ReS pi G(s)/s (2.2) 

i 

where : 
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ReS . = Lim G(s)/ s (s-P.) 

A s+Pi 1 



and P^ are the poles of G(s)/s. 

The function H(z) has generically n-1 zeroes. For particu- 
lar values of the sampling period some zeroes may, however, 
go to infinity, or they may be cancelled by poles, i.e., 
hidden modes. Hidden modes should not be considered as zeroes 
of H(z) . In fact, there are in general no simple closed form 
expressions for the zeroes of H. The limiting cases for small 
or large sampling periods can be characterized. That is ex- 
plained in the following theorem. The major steps in the 
proof are given by [Ref. 9] . 

Theorem 1 : 

Let G(s) be a rational function 



( s — z^ ) (s-z 2 ) . . . (s-z m ) 

(s-p-j^) (s-p 2 ) ••• ( S- P n ) 



(2.3) 



and H(z) the corresponding pulse transfer function. Assume 

that m < n. Then as the sampling period T -* 0 , m zeroes of 

H(z) go to 1 as exp(z^T) and the remaining n-m-1 zeroes of 

H ( z ) go to the zeroes of B (z) where B (z) is the polynomial 

n-m n 

defined as 



B ( z ) 
n 



k v, n_1 

b. n + 

1 z 



H r, 11-2 

b~n + 

2 z 



+ b 



n 



n 



(2.4) 



and 
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n 



(2.5) 



k 



k-1 .n ,n+£ 
l k-£ 



b 



k 



I (-D 



k = 1 



z . * * , n 



1=1 



The following results can be observed from the previous 



1. The limiting zeroes of a pulse transfer function depend 
critically on the pole excess of the corresponding 
continuous time system. 

2. A continuous time system with a pole excess larger 
than two will always give a pulse transfer function 
with zeroes outside the unit disc provided that the 
sampling period is sufficiently short. This may 
happen for quite reasonable sampling periods, and 
sampled data systems with unstable inverses are thus 
quite common. 

An example is given to illustrate the above discussion. 

Example 2.1 : 

Consider a system with transfer function 



The corresponding pulse transfer function can be computed as 



theorem: 



G (s ) 



1 




H (z) 




where : 



b 



1 



1 - (1 +T +T 2 /2)e 



b 



2 



(-2 + T +T 2 /2)e T + (2 +T -T 2 /2)e 
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b 



3 



(1 -T +T 2 /2)e 



e 



-3T 



H(z) has a zero outside the unit disc if 0 < T < 1.8399. 

The result of Theorem 1 can be understood by studying the 
behavior of the continuous time transfer function G(s) = s n 
and its pulse transfer function. The reason why Theorem 1 
holds is because for a sampling period T small every system 
behaves like G(s) = 1/s 11 m with the number of poles n and 
number of zeroes m. 

The pulse transfer function corresponding to G(s) = s n 
is given by 



where B^tz) is given in Equation (2.4). 

The polynomials B n are found for some values of n using 
the program given in Appendix A. 



H ( z ) 



T n B (z) 
n 



( 2 . 6 ) 



n! (z-l) n 



B x (z) 



1 




z +1 




2 

z + 4 z + 1 




z 3 +llz 2 + 1 1 z +1 



B 5 (z) 



z 4 + 2 6z 3 + 66z 2 + 2 6z +1 
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The polynomials B n have zeroes outside or on the unit 
circle for n > 2. The unstable zeroes are given in Table 2.1. 



TABLE 2.1 



UNSTABLE ZEROES OF B (z) 



n 



Unstable Zeroes of B n (z) 



2 



-1 



3 



-3.732 



4 



-1, -9.899 



5 



-2.322, -23.20 



Therefore, it can be noticed that there are continuous 
time systems with stable zeroes such that the corresponding 
pulse transfer function has unstable zeroes. 

It is possible to give a complete characterization of the 
zeroes of the pulse transfer function for small sampling 
periods. A similar result for large sampling periods is 
given by Theorem 2. 

Theorem 2 : 

Let G(s) be a strictly proper rational transfer function 
with G(0) ^ 0 and ReP^ < 0. Then all zeroes of the pulse 
transfer function go to zero as the sampling period T goes 
to infinity. 
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A lower limit on T for stable zeroes was obtained in 
[Ref. 14], here stated only for simple poles. 

T > i £n [2A(n+l) ] (2.7) 



where : 



S = - max ReP . > 0 

l 



( 2 . 8 ) 



and 



A = max I A./ G ( 0 ) | (2.9) 

i 



If G(0) = 0 it follows from Equation (2.10); 



n P i T 

H ( z ) = G(0) z" 1 + (l-z _1 ) l h L - 5 - — • 12,1 

i=l l 

z - 

The corresponding pulse transfer function has one zero at 
z = 1. The behavior of the rest of the zeroes may be more 
complex as is shown by Example 2.2. 

Example 2.2 ; 

Let continuous time transfer function be 



G ( s ) - ^ 

[ (s+ir +1] (s+2) 



Then the corresponding pulse transfer function is 
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